SPECTRAL PROPERTIES OF THE 2D HOLSTEIN t-J MODEL 
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Employing the Lanczos algorithm in combination with a kernel polynomial moment expansion (KPM) 
and the maximum entropy method (MEM), we show a way of calculating charge and spin excitations 
in the Holstein t-J model, including the full quantum nature of phonons. To analyze polaron band 
formation we evaluate the hole spectral function for a wide range of clcctron-phonon coupling strengths. 
For the first time, we present results for the optical conductivity of the 2D Holstein t-J model. 



Polaronic features of dopant -induced charge car- 
riers have been observed in the isostructural 
copper-based and nickel-based charge-transfer 
oxides La 2 - x Sr x [Cu, Ni]04 +y [1]. 

Studying (bi) polaron effects in such strongly 
coupled electron-phonon (EP) systems, the Hol- 
stein t-J model (HtJM) has recently attracted 
much attention [2] . The HtJM Hamiltonian reads 

h = mo Q y, + - ( & I + b i) ~ h i 

i i 



The first two terms take into account the phonon 
part and the EP interaction, respectively, whereas 
the last two terms represent the standard t-J 
model acting in a Hilbert space without double 
occupancy. In (1), doped holes [hi = 1 — fij) 
are coupled locally to a dispersionsless optical 
phonon mode (e p - EP coupling constant, ujq - 
bare phonon frequency). 

In this contribution, we investigate the HtJM 
by performing exact diagonalizations on a square 
ten-site lattice, where the phonon degrees of free- 
dom are treated within a well-controlled Hilbert 
space truncation procedure [3]. To obtain infor- 
mation about dynamical properties of the model 
under consideration, we combine the Lanczos al- 
gorithm with the KPM and MEM approaches [4] . 

In order to address the problem of polaron 
formation in an antiferromagnetic correlated spin 
background, we have calculated the X-resolved 
spectral function Ag{E) for a single dynamical 
hole at J = 0.4 (energies in units of f) [5]. The 



positions of the two lowest peaks of Ag(E), de- 
noted by Eq/x(K), are displayed as a function 



-5.0 



-7.0 



-9.0 



-5.0 



-7.0 



-9.0 



" 0.0 1.0 2.0 3.0 4.0 



□ 



(a) 

h(D =0.8 



O (0,0) 

□ (3,l)it/5 

A (2,4)it/5 

<l (JI.JI) 



0.0 1.0 2.0 3.0 4.0 

e. 



•- 




0,4 




< 


go.2 






% 




0.0 




i 


( 


4 


• 


<] 



(b) 

h(0=3.0 



I I ! 



0.0 1.0 2.0 3.0 4.0 



0.0 1.0 2.0 3.0 4.0 

E 



Fig. 1: Polaron band formation in the 2D HtJM. 
The wave function renormalization factors, 

Zo(K) oc EJ(^ N ' 1 \k)\c^M N \0))\ 2 , are 
given as a function of e p in the insets. 
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Fig. 2: Optical conductivity a^.% 9 (ui) and w-integrated spectral weight in the dissipative part of a re9 , 

= dcu' a™ 9 (to') , for the 2D single-hole HtJM with 12 phonons (periodic boundary conditions), 
of EP coupling strength in Fig. 1 for the al- 



lowed K vectors at Uloq = 0.8 and 3.0. As ex- 
pected, in the very weak EP coupling limit the 
low-lying excitations are t-J hole-quasiparticles 
weakly dressed by phonons. For low and interme- 
diate phonon frequencies, the energy to excite one 
phonon lies inside the quasiparticle band of the 
pure t-J model. Thus at arbitrarily small e p , pre- 
dominantly phononic states with a small admix- 
ture of electronic character enter the low-energy 
spectrum in all i?-sectors [cf. the region about 
E n ~ —5.2 in the inset of (a)]. With increasing 
e p a strong mixing of holes and phonons takes 
place, whereby both quantum objects completely 
loose their own identity, and finally an extremely 
narrow well-separated polaron band is formed at 
large e p . This scenario is corroborated by the be- 
havior of the X-dependent renormalization factor 
Zq(K) shown in the upper insets, which can be 
taken as a measure of the "electronic" contribu- 
tion to the polaronic quasiparticle (see insets) . As 
can be seen from Fig. 1 (b), the phonon induced 
band renormalization is weakened in the non- 
adiabatic regime, where retardation and multi- 
phonon effects arc of minor importantance. 

To discuss the influence of the EP coupling 
on the optical response of the system, let us eval- 
uate the regular part of the optical conductivity 
at finite energy transfer u, 
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Results for g t ^ 9 {lo) are presented in Fig. 2 for 
s p = 0.1 (a), 2 (b), and 4 (c). In the weak 
EP coupling region, we recover the main features 
of the optical conductivity of the t-J model (in- 



set Fig 2. (a)), i.e., (i) an "anomalous" broad 
mid-infrared absorption band [ J < ui < 2t] , sepa- 
rated from the Drude peak [D 5(ui); not shown] 
by a "pseudo-gap" ~ J, and (ii) an "incoherent" 
tail up to u ~ It. At larger e p , we observe a 
redistribution of spectral weight to higher ener- 
gies (cf. ASu), which is much more pronounced 
in the adiabatic regime. In particular, the transi- 
tion to the (hole) polaron state is accompanied by 
the development of a broad maximum in o r ^ 9 (io) 
at lo < 2e pi whereas the optical response becomes 
strongly suppressed at low lo. Most notably, 
g t ^ 9 (lo) has an highly asymmetric lineshape at in- 
termediate frequencies and coupling strengths as 
observed, e.g., for Lao.gSr .iNi04 [1]. This effect 
can be traced back to a rather broad ground- 
state phonon distribution function obtained for 
e p = 2, 4 and hcoo = 0.8 [5]. Contrary, in the anti- 
adiabatic limit, the "electronic" lineshape is much 
less affected, but g t ^ 9 {lo) shows additional super- 
structures corresponding to "interband" transi- 
tions between t-J-like absorption bands with dif- 
ferent number of phonons [see Fig 2 (c), inset]. 
These results clearly demonstrate the complex in- 
terplay of electron and EP correlation effects. 
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